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Abstract 

We consider the phase behavior of polymeric systems by calculating the structure 
factors beyond the Random Phase Approximation. The effect of this correction 
to the mean-field RPA structure factor is shown to be important in the case of 
coulombic systems. Two examples are given: simple electrolytes and mixtures of 
incompatible oppositely charged polyelectrolytes. In this last case, all former studies 
predicted an enhancement of compatibility for increasing charge densities; we also 
describe the complexation transition between the polyelectrolytes. We determine a 
phase diagram of the polyelctrolyte mixture that includes both complexation and 
incompatibility. 

1 Introduction 

Most polymer blends are incompatible and phase separate in a wide range of their phase 
diagram into two phases each containing essentially one of the components of the mixture. 
The physical origin of the segregation of the polymers is their low mixing entropy of order 
where N is the number of monomers per chain, which cannot overcome any weak en- 
thalpic interaction. A slight dilution of the two polymers in a solvent does not change the 
quaUtative picture of the phase separation, except that the two phases in coexistence are 
diluted. The chemical mismatch between two polymers is measured by a Flory-Huggins 
parameter xab that gives the contribution to the free energy Fenthaipic = kT xab <pA(pB 



1 



where the monomer concentrations of polymers A and B are denoted by (pA and 0b . For 
most polymer pairs, xab is positive leading to incompatibility. There exists however some 
cases where xab is negative and where the polymers are compatible and form a single 
homogeneous phase 0. A way to enhance the compatibility in solution is to add some 
electrical charges on the polymers[0]. The simplest case is a mixture between a polyelec- 
trolyte and a neutral polymer. In addition to the ionic charges along the polyelectrolyte 
backbone, the solution contains small counterions that insure the electroneutrality. Their 
entropy is sufficient to overcome the enthalpic interaction between the backbones since it 
is proportional to the monomer concentration independent of the polymer length: a phase 
separation is unfavorable as the counterions loose too much entropy by being confined in 
one of the phases. One can also directly take advantage of the electrostatic interaction 
by charging both polymers with ions of opposite charges. Such experiments have been 
performed by Djadoun et a/. Q by mixing acidic and basic copolymers. If the number 
of charged groups along the copolymers is large enough, a single homogeneous phase is 
observed. The corresponding ternary phase diagram is sketched in Figure l.a. It is qual- 
itatively well described by the theory of Khokhlov et al. for low charge densities along 
the polymer. For higher polymeric charges, the shape of the phase diagram is completely 
different, and is sketched in Figure l.b: the two polymers form complexes leading to a 
polymer-solvent type of phase separation. In this case the Flory parameter is no longer 
relevant since the phase behavior is governed by electrostatics. 

Polyelectrolyte complexes have been studied extensively for thirty years |^ . Nevertheless 
a clear theoretical description is not yet available. A so-called "host" model has been 
proposed by Nordmeier et a/. 0; the complexation is described as an ion binding between 
the polymers. This model is also invoked in a number of experiments 0] and should be 
realistic for highly charged semi-fiexible polyelectrolytes. The first attempt to describe 
complexes between fully fiexible weakly charged polyelectrolytes is due to Borue et al. . 
In this approach, the complex is stabilised by attractive interactions induced by charge 
fiuctuations. We have used this model recently to describe the complexation between 
adjacent layers in polyelectrolyte multilayers 0. 

Our purpose in this paper is to describe with the same model the transition from an in- 
compatibility phase diagram (Fig. la) to a complexation phase diagram (Fig. lb) when the 
charge density along the polymer chains is increased. The general formalism is introduced 
in section 0. It relies on a field theoretical formulation of the Random Phase Approxima- 
tion. In section ^, we apply the method to a simple electrolyte as an example of the new 
results brought by the theory. Then it is applied to a charged polymer mixture where the 
two types of chains show chemical mismatch. We discuss the calculated structure factors 
and the phase diagram of the solution in section 4^. In the last section, we present some 
concluding remarks and discuss some possible issues. 
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2 General theoretical framework 



Our aim in this section is to introduce the formahsm that we use throughout the paper. 
We introduce a Hamiltonian H[(j){r)] that can be written as a functional of the local field 
(below is the local concentration or the composition of the polymer mixture). The 
thermodynamic properties of this system are given by the partition function lir 
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where we introduce an external field h(r). For the rest of this paper we use kT as our 
energy unit {kT = 1) to simplify the notations. The ra-point correlation functions of the 
field are: 

G(r,. ... . . . r„) . . . . ^(r„)) = ^ ,,(,^)f,f,^)'^'',,(,„) \ ,.„^„ (2) 

In this paper we focus on the two-point correlation function (the structure factor). It 
can be calculated as a derivative of the thermodynamic potential r[0(r)] which is the 
Legendre transform of the free energy F[h{r)] = — log Z[h{r)]: 

r[0(r)] = F[h{r)] + J dr/i(r)0(r) (3) 



(4) 



We call 0(r) the thermodynamic average of the field 0(r). The correlation function is 
obtained as: 

G-i(r,,r,) = J^SWL 
50(ri)50(r2 

For any realistic Hamiltonian, the calculation of the path integral involved in Eq.|l] cannot 
be performed exactly. The simplest approximation is to compute the partition function at 
the saddle point of its integrand: this is the mean field approximation. It neglects the con- 
tribution of all other paths to the integral in the partition function. This approximation 
leads to rather poor results when correlation effects are important. One can improve the 
approximation by taking into account the contribution of the fluctuations around the sad- 
dle point path at the gaussian level. This is the so-called Random Phase Approximation 
(RPA). We present here a standard route to derive the RPA thermodynamic potential. 
The field can be decomposed into two parts 0(r) = 0mf(i") + '^0(r) where the mean field 
part stems from the minimisation equation: 

d(f)[r) 4>MF 

Upon expansion of the Hamiltonian around the saddle point at second order in the fluc- 
tuations of the field, the partition function is rewritten as: 

Z[h{v)] = ZMF[h{r)] j D6(/)ir)exp[-^ j dr,dr2 6(Pir,) G],],^{r,,r,) 6<P{r,)] (6) 
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where we introduce the mean field partition function and the mean field two-point corre- 
lation function: 



[Mr) 



const. exp[-i^[0A/i;'(r)] + / (ir/i(r)0MF(r)] 



0AfF{r) 



(7) 



Note that Gj^p^ depends a priori on the mean field path (^mf- One can then calculate 
the RPA thermodynamic potential: 



r[0(r)] = i/[0(r)]-log 



L'50(r)exp[-i j dridr2 5(j){ri) Gpp^{ri,r2) 5(j){r2) 



with ^(r) = -M, 

inverse structure factor reads 

G ^(ri,r2) = GppJ^{rl,r2) 



Using Eq.|, the inverse two-point correlation function, or the 



(50(ri)50(r2) 



log 



D50(r)exp[-i / c/r'dr"(50(r')G^p^(r',r")50(r") 



(9) 



The RPA formalism is of common use in polymer physics when one has to deal with rather 
dense systems where concentration fiuctuations are known to be weak0,[0. One can 
calculate both the thermodynamic potential and the structure factors. The free energy is 
useful to discuss macroscopic phase separation, but the knowledge of the structure factor 
is necessary if one wants to study properties at a finite lenth scale such as a micro-phase 
separation. For this purpose, several theoretical approaches only consider the first term 
of Eq.p, which corresponds to a mean field approximation for the inverse structure factor: 
fluctuation effects are neglected. In most cases, this is sufficient to discuss the phase 
separations of interest. 

However for some coulombic systems, correlation effects might be relevant. Consider for 
example the case of a simple monovalent electrolyte. We show in the next section that 
neglecting fluctuations in the RPA calculation of the structure factor, which is the spirit 
of the Debye-Hiickel theory of electrolytes (DH), leads to a constant density-density struc- 
ture factor that does not satisfy all the sum rules imposed by the electroneutrality. This 
observation is at the origin of the Generalized Debye-Hiickel (GDH) theory of Lee and 
Fisher which produces a g-dependent structure factor | ]12| |. 

One can qualitatively understand why a mean fleld approximation misses some rele- 
vant properties in coulombic system: because of the electroneutrality condition, no pure 
coulombic term can appear in a mean fleld thermodynamic potential in the absence of 
external fleld, since fluctuations around the electroneutral state are not included. 
The RPA calculation of the thermodynamic potential is also called the one loop approxi- 
mation, since it corresponds to a resumation of all the one loop diagrams of the associated 



4 



field theory for tliis potential. In tlie following, we call the result of Eq.^ the one loop 
approximation (OLA) of the inverse structure factor, to distinguish it from the standard 
RPA approximation, that includes only the zeroth order term in a loop expansion. 



3 Structure factor of a simple electrolyte 

We consider in this section a simple monovalent electrolyte solution in water composed of 
discrete point-like ions embedded in a continuous dielectric background of permittivity e. 
For a sake of simplicity, we neglect the effect of the hard core potential since we only want 
to show the effect of the concentration fluctuations and not to discuss quantitatively the 
phase diagram. Simple electrolytes have been investigated quite recently by Netz et al. 
||T3| and Frusawa et al. independently using field theoretical methods. The N-particle 



problem has been mapped exactly onto a field theory with two fields which correspond 
to collective variables: the first one is a density field and the second one is an auxiliary 
field which acts as a complex potential on the ions. At the saddle point of the auxiliary 
field, the partition function is written Z = j Dr?,+ (r)Dr?,_(r) exp[— iJ[r?,+ (r), 'n,_(r)]] with 
the Hamiltonian: 

i?[n+(r), n_(r)] = j (ir|ri_|_(r)(logn+(r) — 1) + ri_(r)(logn_(r) — 1) 

^l_B I ^^, {n+{r) - n_{v)){n+{r') - n_{v')) 



(10) 

2 

The Bjerrum length is defined here by = it measures the strength of electrostatic 
interactions. This hamiltonian can be used as the starting point of a density functional 
theory where the equilibrium density fields result from a balance between entropy and 
electrostatics. At this level, we neglect the fluctuations of the auxiliary field and investi- 
gate the effect of composition fluctuations. We show below that we recover this way the 
results of the Generalized Debye-Hiickel theory if we include the composition fluctuations 
at the level of a one loop approximation, and the classical results of the Debye-Hiickel 
theory if we neglect composition fluctuations. Let no-(r) = ?7.a-(r) + 5n{Y) with a = 
Then using the definition of section ^ the partition function is: 

Z = Zmf j -D5n_|_(r)D5ri_(r) exp[— Aif[?2+, n_, (11) 

where the effective hamiltonian is defined as 



The RPA inverse structure factor matrix is given by: 

S{r-r') , inlg 



47rZf 



'-'RPAa,a'\'^^'^ ) - A-kIb 5(r-r') , inlg ^^'^l 



|r-r'| n_(r) |r-r'| 
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One can easily generalize Eq.^ to a mult i- component system to obtain the inverse structure 
factor matrix in the one loop approximation: 

/ 6^ AH \ / 6AH w 6AH \ i 6AH 6AH \ 

where the average is performed with the gaussian weight exp[— Aif[n_|_, n_, (5ri„]] 
of the fluctuating variables 5n_|_,5n_. The average densities n_|_,n_ are constant in a 
homogeneous phase. Note that the densities must be considered as constant only at the 
step of averages, after the formal derivations with respect to inhomogeneous densities; 
otherwise the contributions of the fluctuations vanishes. 

It is convenient to study the Fourier transform of the structure factors. Each element of 
the matrix is then rewritten as: 

G;,Uq) = G^PA.,.'(q) 



2nlnl, 



i^^Y GrPA a,a' i(ll)GRPA a,a' (q " qi) + J ^^^.^3 Grpa a,a' (qi) 

(15) 



where the gaussian average has been already performed. All the integrals involved in Eq.p!^ 
can be calculated analytically by using Feynman integral techniques. The diagonal terms 
are diverging at high q. This divergence can be removed by substracting the contribution 
of the neutral system (perfect gas). The structure factor matrix of the electrolyte solution 
reads: 

_ /i + ^-AG-^q) -^-AG'-i(q) \ 

G,..(q)- 1^ "_4^_^^-i(q) l+^-AG-Hq) J ^^^^ 

The electroneutrality condition imposes the equality of homogeneous densities = = 
n. The correction term is then given by: 

AG-\q) = ^ arctan ( ^] (17) 



q \2k 

with the Debye-Hiickel length defined by = STclsn- The effect of the fluctuation 
contribution can be drawn from the density-density structure factor that is calculated 
from the partial structure factors Gntot,ntot = Gn+,n+ + Gn_,n- + 2Gn+,n-- 

2n 

Gn,^,,n,M = i_2nAG-\q) ^^^^ 

If we neglect the fluctuations at all orders,the density- density structure factor is constant. 
The fluctuations induce a decrease of the structure factor for increasing wave vectors. 
Notice that the value value at g = is consistent with the Debye-Hiickel polarization 
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energy Fpoi = This is expected since the results of the OLA and the RPA are the 

same for the free energy. 

This natural generalization of the Debye-Hiickel theory to structure factors must be com- 



pared to the Generalised Debye-Hiickel theory proposed by Lee and Fisher ([0. Starting 
from the same constatation on the behaviour of the density-density structure factor in the 
Debye-Hiickel theory, they solve the Debye-Hiickel equation for the electrostatic potential 
with modulated densities. When the amplitude of the modulation vanishes at the end of 
the calculation, they obtain a g-dependent structure factor. Notice that they take into 
account the finite radius of the ions that we neglect here. Nevertheless their procedure is 
very similar to ours. In fact, the correlation length that we can define from the density- 
density structure factor is exactly the same in the limit of vanishing ion radius. More 
precisely the q = susceptibility and the correlation length are in the limit of vanishing 
ion radius: 

x-^(o) = i-^, e = ^ (19) 



in agreement with the results of Lee and Fisher. It is clear from Eqs.|18| and |19| that the 
"One loop approximation" is valid as long as ^ < 1. If it is not the case, one has to 
consider higher orders in the loop expansion. The perturbation parameter of the theory 

is C = 

We do not discuss further the properties of the structure factor because it goes beyond the 
scope of this paper. We note however that simple tests do not reveal any inconsistencies, 
in particular with the Stillinger-Lovett sum rules [O . 



4 Complexation and incompatibility in mixtures of 
oppositely charged polyelectrolytes 

4.1 Inverse structure factor matrix 

In this section we generalise the results obtained in the previous section for simple elec- 
trolytes to polyelectrolyte complexes which are mixtures of polyelectrolytes of opposite 
charges. We use the structure factors to study the stability of the polyelectrolyte mixture 
and compare our results to the experimental work of Djadoun et a/. mentioned in the 
introduction. In this work, the charge of acidic and basic copolymers is progressivley 
increased. The neutral backbones of those copolymers exhibit a strong incompatibility. 
Starting from an incompatibility phase diagram at zero charge, the compatibility range 
(solubilisation) increases with the charge and finally turns into a complexation phase di- 
agram where the two copolymers are in the same neutral dense phase in equilibrium with 
almost pure water because of the strong electrostatic attraction between them. 
In order to describe this system, we consider a mixture of two polyelectrolytes A and B 
having the same linear charge density / but of opposite charges |]TB[. The charge density 
is smeared out along the backbone and the solvent is modelled as a continuous dielectric 
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background of permittivity e. A concentration n+ + n_ = 2n of small point-like ions 
of opposite charges is added to the solution. We study only a symmetric system and 
the monomer concentrations are c+ = c_ = c. Each chain is composed of N monomers 
of size a. If the electrostatic interactions are switched off, the short range interactions 
between polymers are characterised by three Flory parameters Xas,Xbs and Xab- For 
a sake of simplicity, we will only consider both polymers are in a ^ solvent and fix the 
values xas = Xbs = |- The last parameter xab = X characterises the chemical mismatch 
between the two polymers and is varied independently. 

The total Hamiltonian of the solution can then be written as: 

if[c+(r), c_(r), n+(r), n_(r)] = J (ir|ri+(r)(logn_,_(r) — 1) + n_(r)(logn_(r) — 1) 

^a^ /|Vc+(r)|2 |Vc_(r) 



24 V c+(r) c_(r) 



^Ib f^^,P^)p.ir' 



+XC+(r)c_(r)| (20) 

The charge density is defined here as Pzi'r) = /(c+(r) — c_(r)) + n+(r) — n_(r). The 
first term of Eq.^ is the translational entropy of the small ions. The second term is 
the conformational entropy of the polymer chains in the limit of infinite chain length. In 
this limit the conformational entropy can be calculated from a ground state dominance 



approximation and it is given by the Lifshitz formula [17]. In the limit of infinite chain 
length, the translational entropy of the polymer chains is negligibly small. The third term 
is the electrostatic interaction between the charged components of the solution. Finally the 
last term models the enthalpic interactions that give rise to the incompatibility between 
the polymers. 

From this model Hamiltonian we want to calculate various structure factors in order to 
discuss the stability limits of the solution and its phase diagram. As for simple electrolytes, 
the calculation is performed at the level of the one loop approximation presented in section 



For a sake of simplicity, we include the enthalpic interaction between the polymers only 
at the mean field level, or equivalently at the RPA level in the calculation of the structure 
factor. We believe that this approach gives the correct physical result because we expect 
fluctuations to be dominated by the connectivity of the chains. This calculation can also 
be used as a starting point to study various systems involving gaussian polymer chains in 
external potentials (in order to take into account excluded volume or solvent effects for 
example). Note that the contribution of the Flory parameter to the fluctuation can also 
be computed but the calculations are lengthy and much more tedious. 
We first ignore the short range interactions, x = 0- order to find the fluctuation 
contributions to the structure factor, let Co-(r) = Co-(r) +5co-(r) and ^^.(r) = n^{r) + 6n^{r). 
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By expanding the Hamiltonian up to second order in the composition fluctuations, we 
calculate the thermodynamic potential and the resulting inverse structure factor. Within 
the one loop approximation, the inverse structure factor reads: 

/ 52 A/7 \ . / 5^H \ I 5^H \ I 5^H 5 AH 



^\5p,{v)5-p,. (rO / + \ 5p.(r) / \ dp,. (r') / \ 5p,{r) bp^, (r') / ^^^^ 

with the short-hand notation p^, = (c+, C-,n+,n^). The quantity AH is the second order 
term in the expansion of the Hamiltonian H. As in the previous section, the average is 
performed with the gaussian weight exp[— Aif] of the fluctuating densities. The average 
densities p can in the calculation of the averages be set equal to their uniform values in 
the solution, c and n respectively for the monomers and the small ions. The RPA inverse 
structure factor has been already calculated by several authors , it is given in reciprocal 
space by: 



I TOT^ ' q2 „2 „2 „2 \ 

(2 " " 



^RPAp, o-.cr' 



12c + f 



i-KlBf _MbI 

9 9 - ~r 9 9 

^ gZ gA 

\ (p (p n q'^ / 



(22) 



We give below the values of the corrections to the RPA in Eq.^ for the various components 
of the inverse structure factor matrix. 
The monomer-monomer components read: 

T^[-^-ir [ JTT^ i^l-i^- ^l) - Grp Aca,a'{q.l)G RPAca,a' (ci- q.l) 

576 c^c;, J {2TTy 
, 48 f d^qi 



(23) 



The correction term to the RPA for the salt-salt components of the inverse structure 
factor matrix are given by the same formula Eq.|l^ as for a simple electrolyte: 

^g;X.' = 

ni J {^TlY 

(24) 



2nlnl, j (27r)3 



Note however that the RPA structure factor is different from the simple electrolyte case. 
Finally the corrections to the crossed salt-monomer terms of the inverse structure factor 
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matrix are given by: 

= 24c%2, J (2^ ^"^'-^"^ ~ "^'^ ~ "l'^ GRPAc„,nA^l)GRPAc.^nA^ " ^l) (25) 

All the integrals involved in Eqs.^, ^ and ^ can be calculated using Feynman integrals 
techniques. The analytical expressions of these integrals are reported in appendix A. 

We can therefrom calculate the inverse structure factor at the level of the one loop 
approximation. As the inverse structure factor matrix is symmetric, we only give the 
independent terms: 

2 2 A J f2 f2 

G-oIac,,c, = GoIacc. = \^ + ^^ + -^^D-f"f^^'D (26) 
GoIa.,,. = -^-fm (27) 

GoLAn+,n+ = GoLAn-,n- = ^^^) (28) 

GoLn,,n_ = -^-^(^) (29) 

GoLc,,n, = GoLc.,n_ = ^ " f 7(g) (30) 

GoLc„n_ = GoLc_,n, = ' fH^) (31) 

The quantities a, /3, 7, 5 are related to the Feynman integrals of appendix A by: 

"(g) = — —h{q) 

m = (32) 



4.2 Phase diagram and discussion of the experimental results 

We have obtained in the preceding section the inverse structure factor matrix in the 
particular basis {c+, c_, n+, In order to discuss the phase diagram of the system, it is 
more convenient to calculate the structure factor matrix for the variables {ctot, Cz, ntot, n^} 
denoting respectively the total monomer density {ctot = c+ + c_), the total monomer 
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charge density (c^ = /(c+ — c_)), the total small ion density {ritot = n+ + n_) and the 
total small ion charge density {uz = n+ — n^). Indeed the density- density structure 
factor is essential to predict a phase separation such a complexation which is a polymer- 
solvent phase separation; it is associated to a singular behavior of this structure factor. 
On the other hand, the incompatibility between the chains, is a polymer-polymer phase 
separation, and is associated to a singularity in the charge-charge structure factor, since 
the polymer charge is the only way to distinguish the two types of chains within our 
model. Previous studies have focused on this last structure factor because it exhibited a 
singularity while the density- density structure factor was regular. We show below that 
the effect of the composition fluctuations is to induce an attraction between the chains 
that counterbalances in some parameter range the chemical mismatch between the chains. 
This effect was ignored in all the previous studies that stay at the RPA or mean field level 
for the structure factor. 

The change of variables and the inversion of the matrix Gp^^j^^^q) can be performed 
analytically. At this point we add at the RPA level the Flory parameter x- The structure 
factor bloc matrix for the polymer variables then reads: 

GoLAc^,c^,{l) = 

^ ^ — \ 

' (33) 



^(4#+/'°(9)-2/*c/3(9)) + 



n i 

where a = {tot,z). In this last equation, we have neglected higher order contributions 
obtained in the matrix inversion involving the square of correction terms since these orders 
in the loop expansion have been neglected from the very beginning. 
If we neglect a{q) and /5(g) in the structure factor, we recover the results of Khokhlov 
et their discussion is limited to the charge-charge structure factor. We present 

here a similar analysis by omitting first the effect of the density-density correlations on 
the phase diagram. The phase diagram can be drawn in a very transparent and elegant 
way by using the dimensionless variables introduced first by Borue et al^. So far, we 
have introduced two characteristic lengths: the Debye-Hiickel screening length and 
the correlation length of a dense polyelectrolyte mixture (g^ = 2^ilB]_cy r^-^^^ gj,g^ 
length is associated to the screening of electrostatic interactions by the small ions and 
the second length is associated to the screening of the same interactions by the polymers. 
The enthalpic interactions between the neutral backbones of the chains are measured by 
the Flory parameter x to which we can associate the mean field length C,x defined by : 

- §^ (34) 

The physical behaviour of the system is then governed by the following dimensionless 
ratios: 
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'f. 



t = ^ (35) 



q* 

The choice of as the unit length is motivated by the fact that we are expecting the 
connectivities of the chains to play a dominant role. 

The stability of the system against the polymer-polymer phase separation is ensured if 
1/GczCz(q') is positive for all values of q. The spinodal line giving the stability limit is 
obtained when this condition is violated. If the first mode for which 1 / Gc^c^ (?) becomes 
negative is go 7^ 0, the solution shows a microphase separation. On the contrary, for go = 
there is macrophase separation. The partial phase diagram is shown in Figure 2. It does 
not take into account yet instabilities of the density-density structure factor. An analysis 
based on the stability of the homogeneous phase can only give the equation of the spinodal 
lines but it cannot be used to predict the precise morphology of the different phases after 
the transition. A complete phase diagram can only be obtained by calculating the free 
energy of the various possible phases and by looking for the phases that minimise it. In 
this paper wc study only the stability limits of the homogeneous phase. 
We checked by numerical evaluation of the structure factor with the exact formula that 
the correction term a{q) can be neglected. If s > 1 the system is homogeneous for st < 1, 
otherwise it is macroscopically phase separated: chains of opposite charge are segregated 
into two dilute phases. This is due to the incompatibility of the polymers. If s < 1 it is 
possible to have microphase separation. The spinodal line for this transition is given by 
t = 2 — s. The wave vector at the transition is q^ = ql^- It is related to the wavelength 
of modulation of the mesophase and it tends to infinity as we approach the macroscopic 
spinodal line st — 1. Notice that in the limit of high ionic strength (s ^ 1) which is the 
relevant case for most experiments, no microphase separation is expected. 

The full phase diagram includes possible instabilities associated to the density-density 
structure factor. The numerical analysis of this structure factor shows that the correction 
term a{q) can again be neglected and that the term /3(g) can be approximated by a 
constant /3{q) fti /3(g = 0). The result of this analysis is that the composition fluctuations 
induce a macroscopic complexation between the oppositely charged polyelectrolytes when: 

where we introduce the following dimensionless ratios: 

^mesh /o^7^ 

u = (37) 
q* 

3 

^mesh — 75 5~ (38) 
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The length C,mesh is the mesh size of a gaussian transient network (gaussian chains in 
semi-dilute regime). This length is relevant to discuss the validity limits of the OLA or 
the RPA calculations. The main assumption of those approximations is that the chains 
are gaussian at all length scales. For simple polyelectrolyte solutions, this remains true 
as long as the electrostatic blob size Q is larger than the mesh size of the semi-dilute 
solution because the electrostatic interactions are screened by the transient network at 
large length scales|]l6|. Up to irrelevant prefactors this condition is equivalent to u < 1. 

;4/3 2/3 

In terms of polymer concentration this threshold is written ca ^ ^ a-^r^ ■ Therefore the 
approximation used in this paper is valid at moderately high concentration. 
The phase diagram of the system is presented in Figure 3. Note that there is no microphase 
separation associated to the polyelectrolyte complexation. Although the variables used 
to describe the phase diagram look very physical, they are not well suited to describe the 
different phases in terms of concentration. In Figure 4 we translate the t — s diagram 
into a n — c diagram (all the other parameters being kept constant). The equations of the 
various spinodal lines are reported in Appendix B. 

If we now come back to the experiments of Djadoun et al. discussed above, our model 
is able to describe qualitatively the transition from incompatibility to compatibility and 
finally complexation for increasing polymeric charges. Indeed the parameters t and s scale 
hke: 

^ - hr-, (39) 

so that the two parameters have the same dependence in the polymer charge density, 
except for prefactors. If we increase /, we move in the phase diagram on a line represented 
in Figure 3. Starting at low /, the solution is in the incompatibility range; increasing 
/ pushes it into the compatibility range and finally complexation takes place. A very 
rough model of the complex can be made by describing it as a collapsed polymer |0. 
This argument also shows that incompatibility becomes an irrelevant parameter once the 
complex is formed. 

In real systems, the polymers are not symmetric in length and charge, but we believe 
that this will not affect significantly the general shape of our proposed phase diagram: 
the neutral dense phases should be similar in composition but we expect the supernatant 
to be composed of charged soluble species. The description of the supernatant will be the 
topic of a forthcoming paper. 



^The electrostatic blob size is the length scale at which the typical electrostatic energy of a subunit 
composed of {^^f^ becomes larger than the thermal energy. More precisely kT w ^'^^-^ ^ (^) 
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5 Concluding remarks 



We have proposed in this article description of polymer mixtures in a solvent that goes 
beyond the classical Random Phase Approximation and that allows to describe simultane- 
ously in a certain range of parameters both the incompatibility between the two polymers 
and the formation of polyelectrolyte complexes when the polymers are charged. This one 
loop approximation has been used to calculate the structure factors. The Random Phase 
Approximation is commonly used to discuss phase stability in polymer systems but in the 
calculation of the structure factors, it misses important correlation effects. As pointed 



out by Haronska et al. |]T8| , the RPA is able to predict only repulsive electrostatic contri- 
butions to the mean field effective potential between polymers. Including fluctuations at 
the order of one loop allows for attractive effective interactions. Using a One Loop Ap- 
proximation at the level of the structure factors, we were able to include the effect of the 
fluctuations for two systems: simple electrolytes and mixtures of incompatible oppositely 
charged polyelectrolytes. For simple electrolytes, we obtain a natural generalization of 
the Debye-Hiickel theory which leads to results similar to the one proposed by Lee and 
Fisher. For polyelectrolytes blends, we were able to describe within a single model the ex- 
periments showing a transition from an incompatibility phase diagram to a complexation 
phase diagram. The transition towards polyelectrolyte complexation is predicted to take 
place in the range of validity of the OLA, namely in the limit of weakly charged chains 
The main limitation of this work is that for high enough charge densities the one loop 
approximation breaks down: the system becomes strongly fluctuating and one should use 
another formalism. A model based on the formation of ion pairs seems appropriate in 
this case; whe hope to be able to set up such a model in a future publication. 

Acknowledgments: We are grateful for discussions to I. Erukhimovitch (Moscow 
state University). This research was supported by the Deutsche Forschung Gemeinschaft 
through the Schwerpunkt program " Polyelektrolyte" . 

Appendix A 

A.l Exact expression for the monomer-monomer correction terms 

By replacing the expression of the RPA structure factor into Eq.|23|, we obtain for the 
off-diagonal term: 

AG--.-l(l?^)%.(.) ,41) 
/i(q) is a Feynamn-like integral defined by: 



/i(q) 



f^^qi (qi-(q-qi) -q^)^ 



^42) 
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with the characteristic length defined by qf = ^^ttZb/ interpreted as the 

correlation length of a dense mixture of oppositely charged polyelectrolytes in the absence 
of small ions . By reducing the integrand of Eq.|4^ to simple fractions we can expressed 
the whole integral as a sum of Feynman diagrams. We give here only the expansion at low 
q because the discussion on the possibility of macro- or microphase separation is relevant 
only in this range, while the exact formula is rather lengthy: 

STcq+q.iq^ + q.)^^'^ 9QT^qlq^{q+ + q_f ^ > 

where the two wave vectors g+, g„ verify q\q^ = q* and (g+ + = + 2g^. 

The diagonal terms exhibit a divergence at high q which can be removed by substract- 
ing the contribution of the neutral system as in section |^. The convergent term is given 
by: 

2At^IbP ^ , , 1 f^SnlBp 

with the new integral: 



^Gc.,.= o l2ici)--( ] /i(q) (44) 



d\i 1 r (qi-(q-qi) -q^)^ (q^ + q?)(q- qi)^ - (q-qi)^ 



(27r)3 q?(q2 - q?) i (q - q02((q - qi)^ + t,^) + qt q2(q2 + + qf 



The low q expansion of this last integral is given by: 



(45) 



6vrg+g-(g+ + g-) 

A. 2 Exact expression for the salt-salt correction terms 

The diagonal and off-diagonal correction terms in the salt-salt matrix are equal after 
substraction of the unphysical divergence and the common value reads: 



AG-l, = -Sn%h{q) (47) 



with the integral: 



T (n]= I i!SL qi(q-qi)' 

^^^^ J (27r)3[q2(q2 + «:2)+^4][(q_q^)2((q_q^)2 + ^2)+^4] I ) 

Although we calculated exactly this integral we do not give here the result because it is 
lengthy and we do not need it for the discussion of the structure factor as it will appear 



clearly in the section |4^ . 
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A. 3 Exact expression for the salt-monomer correction terms 

As in the salt-salt matrix, the diagonal and off-diagonal correction terms of the salt- 
monomer matrix are equal to: 



AG-,,. ^ -'^Ui.) (49) 



with the following integral 



T ( ] = f ^^^1 (q^ -qi.(q- qi)) 

^^^^ J (27r)3[q?(qf + «:2) + ^4][(q_q^)2((q_q^)2 + ^2) + ^4] I ) 

As in the preceding subsection, we will not give the precise evaluation of the integral for 
the sake of clarity. 

Appendix B 

Figure 4 represents the n — c translation of the t — s phase diagram found in the main text. 
We have redefined the small ion concentration by substracting the polymer counterions; 
it is thus such that at = there are still some small ions, the respective counterions 
of the polycations and polyanions. With this convention, the macroscopic spinodal line 
given by = 1 is written as: 

n = ^ - /g (51) 

with the following rescaled parameters: 

h = ho^Ib (52) 
c = co^Ib (53) 

X = 4 (54) 
The spinodal corresponding to microphase separation {t = 2 — s) is given by: 

^ = —aB / + c (55) 



7rV2 y 27cIb 
Finally the spinodal of complexation (t = ^^_^_^)3/2 ) is given by the formula: 

. ^ ^ifP _ _ .... 

At fixed c, n and x, the critical value of the charge density for polyelectrolyte complexation 
is given by the solution of the quadratic Eq.KB| in /. 



16 



References 

[1] p. G. de Gennes Scaling concepts in polymer physics; Cornell University, Ithaca, 
N.Y., 1979. 

[2] T. Kwei, T. Nishi, R. Roberts . Macromolecules, 7, p667, 1974. 

[3] A. R. Khokhlov, I. A. Nyrkova. Macromolecules, 25, pl493, 1992. 

[4] S. Djadoun, R. N. Goldberg H. Morawetz . Macromolecules, 10, pl015, 1977. 

[5] See for example the references inside B.Philipp, H. Dautzenberg, K. J. Linow, J. 
Kotz, W. Dawydoff . Prog. Polym. Sci. , 14, p823, 1998. 

[6] E. Nordmeier, R Beyer . J. Polym. Sci. Part B, 37, p335, 1999 ; T. Schindler, E. 
Nordmeier Macromol. Chem. Phys., 198, pl943, 1997. 

[7] E. Tsuchida . J.M.S.-Pure Appl. Chem. , A31, pi , 1994. 

[8] V. Y. Borue, I. Y. Erukhimovich . Macromolecules 21, p3240,1988 ; V. Y. Borue, I. 
Y. Erukhimovich . Macromolecules 23, p3625,1990. 

[9] M. Castelnovo J. F. Joanny . Langmuir, 16, p7524, 2000. 

[10] P. M. Chaikin, T. C. Lubensky . Principles of condensed matter physics. Cambridge 
University Press, Cambridge. 

[11] L. Leibler . Macromolecules 13, pl602, 1980. 

[12] B. P. Lee, M. E. Fisher . Phys. Rev. Lett. 76, p2906, 1996 ; B. P. Lee, M. E. Fisher 
. Europhys. Lett. 39, p611, 1997. 

[13] R. R. Netz, H. Orland . Europhys. Lett. 45, p726, 1999 ; Preprint of proceedings of 
Les Houches school 2000 .'Electrostatic effects in soft matter and in biophysics. 

[14] H. Frusawa, R. Hayakawa . Phys. Rev. E 61,R6079, 2000. 

[15] F. H. Stilhngcr, R. Lovett . J. Chem. Phys. 48, p3858, 1968 ; F. H. Stillinger, R. 
Lovett , J. Chem. Phys. 49, pl991, 1968 . 

[16] J.L. Barrat, J.F. Joanny . Adv. Chem. Phys, 94, pi, 1995. 

[17] J. des Cloizeaux, G. Jannink . Les polymeres en solution. Editions de physique, Les 
Ulis, FRANCE. 

[18] P. Haronska, T. A. Vilgis . Phys. Rev. E 50, p325, 1994. 



17 



Figure Caption 



Figure 1 la.Typical phase diagram for a solution of oppositely charged polyelectrolytes 
at very low charge densities. In a wide range of this phase diagram the polymers 
undergo a polymer-polymer phase separation, the neutrality being ensured by the 
small ions. The thick lines are the tic lines indicating coexisting phases, 
lb. Typical phase diagram for a solution of oppositely charged polyelectrolytes with 
higher charge densities. In a wide range of this phase diagram there is coexistence of 
dense phases containing both polymers and dilute phases with almost no polymers: 
there is thus complexation between the polyelectrolytes. The thick lines are the tie 
lines. 

Figure 2 Partial t—s phase diagram: the different phases are indicated with their relative 
boundaries. This phase diagram can be associated to the representation used in 
Figure la: the polymers are either incompatible or compatible (soluble). 

Figure 3 Complete t — s phase diagram: the effects of composition fluctuations are taken 
into account, resulting in the possibility of complexation for low x- When the charge 
density / is increased, the other parameters being imposed, the system follows 
the double arrow sketched on the phase diagram, from the incompatibility region, 
through the compatibility region and finally through the complexation region. 

Figure 4 Complexation phase diagram in the concentration variables. The value of x is 
fixed (x = 0.001) and the complexation line is shown for various values of / (from 
bottom to top / = 1/50, 1/20, 1/10). 
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